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Abstract 

Bifurcation properties, stability behavior, dynamics, and the heat transfer of convection struc- 
tures in a horizontal fluid layer that is driven away from thermal equilibrium by imposing a vertical 
temperature difference are compared with those resulting from imposing a heat current. In par- 
ticular oscillatory convection that occurs in binary fluid mixtures in the form of travelling and 
standing waves is determined numerically for the two different driving mechanisms. Conditions are 
elucidated under which current driven convection is stable while temperature driven convection is 
unstable. 

PACS numbers: 47.20.-k, 47.20.Ky, 47.54.+r, 47.27.Te 
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Many nonlinear dissipative systems that are driven away from thermal equilibrium show 
selforganization out of an unstructured state: A structured one appears when the driving 
exceeds a critical threshold The driving might be realized by imposing a field gradient 
across the system — e.g., a voltage difference across a semiconductor [2| or a liquid crystal 3j, 
a temperature difference across a fluid layer 4], or a concentration difference in a physical 5|, 
chemical |6], or biological Ml system — which drives a current. Or, alternatively, a current 

nnn 

might be injected at one side of the system ja, |9|, llfj . Now the question is, whether and how 
the dissipative structures that form in response to these two different driving mechanisms 
are related to each other concerning their dynamics, their structure, their stability behavior, 
and their bifurcation properties. 

lv for the case of convection in a horizontal 



Unlike one-component fluids like pure 



We have investigated this question numerica 
layer of a binary fluid like, e.g., ethanol- water 
water, this system shows a surprisingly rich variety of different convection structures already 
at small driving There are spatially extended states of stationary convection 

rolls and of temporally oscillating roll patterns in the form of traveling waves (TWs) and 
of standing waves (SWs) that bifurcate out of the quiescent fluid state. In addition there 
are also spatially localized traveling wave states that compete with extended convection 
structures. 

Here, we focus on convection in the form of straight parallel rolls as they occur, e.g., in 
narrow channels with roll axes perpendicular to the long side walls. We have solved the 
appropriate hydrodynamic field equations Q| with a finite-differences method 3| in a 
vertical x — z cross section through the rolls perpendicular to their axes thus ignoring effects 
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that come from field variations along the roll axes 

Calculations were done for ethanol-water parameters, Lewis number L=0.01 and Prandtl 
number <r=10. Results are presented here for two different separation ratios ^=-0.03 and 
■^=-0.1 [3| for which TW and SW solutions bifurcate sub critically out of the quiescent fluid 
state via a common Hopf bifurcation. Our findings concerning the above posed questions 
are representative also for TWs and SWs at stronger, i.e., more negative Soret coupling 
strength ip. 

The horizontal boundaries at top (z=l) and bottom (z=0) are no-slip, impermeable, 
and perfectly heat conducting, thus enforcing the absence of lateral temperature gradients 
there. Sapphire or copper plates provide a good experimental approximation. Two different 
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experimentally realizable horizontal boundary conditions (be) for the temperature are ex- 
plored here: (i) Dirichlet be of fixed temperatures (constant in space and time) at z=0 and 
z=l with a difference of AT between them and (ii) von Neumann be of fixed total vertical 
heat current at z=0 and Dirichlet be of fixed temperature at z—\. At the impermeable 
boundaries the vertical concentration current vanishes and consequently the local vertical 
heat current reduces there to —d z T Q|. Note that we impose in case (ii) the horizontal 
mean 



Q = -d z T\ z=0 (1) 

of the heat current at the lower side of the fluid layer so that the total heat current injected 
into it is a constant. We shall identify the driving conditions of case (i) by TT for short and 
those of case (ii) by QT. 

Laterally we impose for all fields periodic be with wavelength A=2. This is roughly the 
critical one for onset of oscillatory convection. Moreover, it is often seen also in nonlinear 
convection with rolls of about circular shape. Finally, to determine SW solutions that 
are unstable against horizontal mirror symmetry breaking phase propagation we enforce 
horizontal mirror symmetry, say, at x=0 thereby fixing the phase 21| . 

As control parameter measuring the strength of the driving we use in the TT case the 
relative deviation 

e = AT/AT C - 1 (2) 

from the critical temperature difference AT C for onset of convection. The driving in case 
QT is measured by the relative deviation 

p = Q/Q c - 1 (3) 

from the critical imposed heat current. 

We shall discuss first TW convection and then SW solutions. In Figs. Q and El we show 
for two different ip the bifurcation diagrams of nonlinear relaxed TW states: (i) current Q as 
a function of e for TT driving and (ii) temperature difference AT versus p for the QT case. 
Note that Q as well as AT are constant for TWs. The diagonal line shows the linear diffusive 
relation Q con dlQc = AT conrf /AT c of the quiescent conductive state. It loses stability via a 
Hopf bifurcation at AT C or Q c , respectively. After increasing the driving slightly beyond 
this threshold transient growth of oscillatory convection occurs with increasing Q for the TT 
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case. For QT conditions AT decreases since convection cools the lower boundary. Initially, 
the oscillations have the Hopf frequency. But finally, the TT transient ends in a stationary 
convection state since the TW branch terminates with zero frequency in a stationary solution 



branch already below e=0 for our ^'s |22j |. On the other hand, the QT growth transient 
ends in a relaxed nonlinear TW (lower part of Figs. Q and EJ) . 

The curves of Q/Q c versus e and of AT/AT C versus p in Figs. ^ and El are reflections of 
each other at the diagonal, bisecting conduction line. Note, however, that the transients and 
the stability ranges of the relaxed TWs are different. Concerning the latter, for example, 
the hysteresis interval in p for QT is significantly smaller than the one in e for TT since a 
large portion of unstable TT generated TWs below onset gets stabilized under QT driving. 

In Fig. El we show for the TWs of Fig. ^ bifurcation diagrams of Nusselt number N, 
reduced frequency uj/lu c , and squared maximal vertical velocity u»^ aa , versus the respective 
control parameters. The Nusselt number 

N = Q/Q C ond = (Q/Q C )AT C /AT (4) 

provides the relation 

p={l + e)N-l (5) 

between equivalent control parameters e and p corresponding to reflection at the conductive 
diagonal in Figs. Q and |2 TWs that are generated by TT or QT driving at e- and p- 
values related by J5J) have the same spatiotemporal properties, e.g., the same N,u,w max as 
indicated by the symbols in Fig. El Their stability, however, might differ. 
Eq. (jSJ yields also the relation 

d p A = d e A/[N+(l + e)d e N] (6) 

between the slopes d p A(p) and d t A(e) in the QT and TT bifurcation diagrams of any order 
parameter A (say, N^uj^w^^etc) versus p or e, respectively. Hence, the QT bifurcation 
becomes already tricritical, i.e., it changes from backwards to forwards when the initial slope 
s = d e N(e = 0) of the TT Nusselt number increases beyond —1. In other words, all TT 
driven backwards bifurcating unstable TWs for which s > — 1 can be stabilized by switching 
over to QT driving. 

Note that the relations (j2J) and © hold also for any stationary convection solution so that 
the bifurcation diagrams of Q/Q c versus e and of AT '/ AT C versus p are reflections of each 
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other. Thus, the stabilization effect of QT driving holds also for any stationary convection 
that bifurcates backwards with TT j^J. The (stability) properties of forward bifurcating 
stationary convection remain unchanged when using QT instead of TT conditions. 

In the remainder of this letter we dicuss SW convection. Under TT (QT) dri ving the 
heat current Q (temperature difference AT) oscillates with twice the SW frequency 25]. So, 
in Figs, ^and El we show bifurcation diagrams of the time averages (Q) /Q c and (AT) /AT C , 
respectively. Like for TWs, QT conditions have a stabilizing effect also on SWs. Note, 
however, that the two SW solution branches in these Figs, are not reflections of each other 
at the conduction diagonal. Their spatiotemporal properties differ and the relation (p) = 
(1 + e) (N) — 1 provides only an approximate equivalence between the bifurcation diagrams 
of the order parameters in Fig. H] and El For example, the SWs marked by symbols in Fig. 03 
have the same frequency. But w max differs slightly and so does (N) — i.e., (Q/Q c ) AT C /AT 
for TT driving in comparison with (Q/Q c ) (AT C /AT) for QT driving. Also the oscillations 
of the flow differ slightly [Fig. EJc)]. On the other hand, the oscillation profile of Q(t) differs 
significantly from the one of AT(t) [Fig. Efa)] and also the profile of NTr{t) differs from the 
one of N QT (t) [Fig.Efb)]. 

In summary: Driving convection with a fixed heat current can stabilize SW, TW, and 
stationary states that bifurcate backwards and that are unstable with imposed temperature 
difference. However, irrespective of their stability relaxed TWs for the two driving mecha- 
nisms are simply related to each other. The same holds for stationary solutions. But the 
time evolution of, say, growth transients differ in general. Also SW oscillations driven by a 
constant field gradient differ from those resulting from constant current driving. It would be 
interesting to see how far these different conditions influence the spatio temporal properties 
of convection structures with more complex dynamics. 
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2l|: 5T(x,z,t) = -5T(x + X/2,l-z,t) holds 



for the temperature deviation from the mean. For QT conditions the temperature T(z = 0, t) 
oscillates almost in antiphase to the current at z = 1. 
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FIG. 1: (color online) Bifurcation diagrams of heat current and temperature difference for oscil- 
latory convection at ip = —0.03 subject to different be. The upper left (lower right) part shows 
Q/Qc (AT/AT C ) for TT (QT) driving versus e (p) on the upper (lower) abscissa. For SWs the time 
averages (Q) /Q c and (AT) /AT C , respectively, are plotted. The bisecting line marks the quiescent 
conductive state. Full (dashed) lines and filled (open) symbols denote stable (unstable) states. SW 
solutions were obtained with phase pinning conditions; otherwise they are completely unstable. 
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FIG. 2: (color online) Bifurcation diagrams of heat current and temperature difference as in Fig. ^ 
Here, however, for ip = —0.1. 
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FIG. 3: (color online) TW bifurcation diagrams: (a) Nusselt number N, (b) reduced oscillation 
frequency u/u c , and (c) squared maximal vertical velocity w^ux f° r ^he TWs of Fig. ^ with ip = 
—0.03. The left and right curves refer to TT and QT driving, respectively. Full (dashed) lines 
denote stable (unstable) TWs. Symbols identify examples of equivalent TWs for equivalent control 
parameters e and p. 



February 8, 2008 

1.008 

A 
V 

1.004 h 



n 



3 



0.75 

o 

0.5 

0.25 
1.5 



CM. B 



0.5- 



1 ' 1 

\^ — ^ / 

\ / 
\ / 
\ / 
\ / 
\ / 

■ A ■ ■ 


(a) 




/ X 

/ \ 

/ \ 
/ \ 
/ \ 

^^^^ \ 
^^^^^ 


— ■ 1 ' 

(b) 

- 


: ' 1 ' 1 


— ■ — i — ' . 

(c). 




\ 7 — — - 




\ / 


- 


\ / 




\ / 




\ / 




\ / 




V 


1 



-0.005 

8 = AT/ AT -1 



0.005 0.01 

P = Q/Q -1 



FIG. 4: (color online) SW bifurcation diagrams: (a) time averaged Nusselt number (N), (b) reduced 
oscillation frequency u>/lo c , and (c) squared maximal vertical velocity hj^ for the SWs of Fig. ^ 
with ij) = —0.03. The left and right curves refer to TT and QT driving, respectively. Full (dashed) 
lines denote stable (unstable) SWs. 
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FIG. 5: (color online) SW bifurcation diagrams as in Fig.|l]but for the SWs of Fig.[2]with ifj 
Symbols identify examples of SWs with the same frequency. 
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FIG. 6: (color online) Oscillation profiles of the SWs marked in Fig. [5] by symbols, (a) Current 
and temperature oscillation subject to TT and QT driving, respectively, (b) Nusselt numbers 
N TT (t) = [Q(t)/Q C ]AT C /AT = Q(t)/Q cond and N QT {t) = (Q/Q c )[AT c /AT(t)] = AT cond /AT(t). 
(c) Vertical velocity w at midheight between two rolls. 



